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ABSTRACT 

A  computer  program  is  presented  which  solves  the  search 
planning  problem  for  survivors  at  sea.  The  program  is 
designed  to  utilize  weather  data  as  compiled  by  the  United 
States  Navy  at  its  Fleet  Numerical  Weather  Central, 
Monterey,  California. 

An  investigation  is  also  made  into  the  search  criteria 
used  by  the  United  States  Coast  Guard  in  its  planning 
procedures.  Guidelines  are  given  for  the  use  of  the  square 
search  and  the  Sector  search. 
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I.   INTRODUCTION 

A.   THE  COAST  GUARD  AND  SURVIVOR  SEARCH 

Survivor  search  is  defined  as  the  process  of  planning 
and  executing  a  search  for  survivors  at  sea,  and  the  Coast 
Guard  usually  accomplishes  this  in  the  following  steps: 

1.  Initial  Notification 

This  stage  involves  gathering  all  available  data 
about  persons  or  craft  in  distress,  including 
characteristics  of  the  survivor  craft,  number  of  survivors, 
etc. 

2.  Initial  Distress  Position 

In  this  stage,  the  search  planner  puts  together  all  the 
available  data  concerning  the  survivor  and  makes  a  judgement 
as  to  what  the  most  likely  position  of  the  survivor  was  at 
the  time  the  craft  seased  operating  normally. 

3.  Drift  Calculations 

The  survivor's  suspected  motion  is  plotted  assuming 
drift  in  the  open  ocean,  and  a  position  is  calculated  for 
some  future  time,  usually  the  planned  beginning  of  search 
efforts. 

U.   Calculation  of  the  Area  to  be  Searched 

An  estimate  is  made  of  the  accuracy  of  the 
survivor's  position  as  calculated  in  step  3  and  the  accuracy 
of   navigation  of   the  search   craft.   An   area,  called   the 


Search  Area,     Is  defined  based  on  these  errors. 

5.  Allocation  of  Search  Units 

A  determination  is  made  of  the  resources  (ships, 
planes,  boats/  etc.)  to  be  used  in  the  search,  and  the  way 
in  which  these  resources  will  be  used. 

6.  Execution  of  the  Search 

Searching  is  accomplished  according  to  several 
established  search  patterns,  most  of  which  are  different 
arrangements  of  a  series  of  parallel  paths  within  the  Search 
Area. 

7.  Termination  of  the  Search 

Either  the  survivors  are  located  or,  after  a  certain 
amount  of  searching  is  completed  without  success,  the  search 
is  terminated.  This  stage  includes  reports,  and  the 
evaluation  of  search  effort. 

B.   DOCTRINE 

The  guidance  for  all  these  activities  Is  contained  in 
Coast  Guard  publication  CG-308,  the  National  Search  and 
Rescue  Manual,  (Ref.  1.).  Chapters  3  through  8  of  this 
manual  contain  the  Instructions  for  planning  and  conducting 
a  search  for  survivors  at  sea. 


C.   PURPOSE 

This  thesis  accomplishes  two  things: 

1.  Provides  a  computerized  solution  to  step  3  of  the 
search  planning  process,  calculating  survivor  craft  drift, 
and 

2.  Investigates  the  criteria  used  for  allocation  of 
resources  to  a  search,  steps  k,    5,  and  6. 

It  is  Intended  that  the  results  be  understood  by  the 
search  planner,  so  that  they  may  be  of  some  use  to  him. 


I  I .  THE  SEARCH  PLANNING  PROGRAM 


A.  OBJECTIVE 

The  major  objective  in  designing  a  computer  program  for 
search  planning  was  to  have  a  fast,  accurate  means  of 
obtaining  an  estimate  of  a  survivor's  oceanic  position  at 
the  time  of  the  beginning  of  a  search. 

B.  THE  MANUAL  METHOD 

In  order   to   understand   the    program,   it   is   first 
necessary  to  understand  how  drift  plotting  is  done  manually. 

l.  Example 

A  fishing  vessel  Is  believed  to  have  sunk  at 
latitude  U0-U2.0N,  and  longitude  52-13. OW,  at  0900z  on  22 
June.  The  probable  error  of  this  position  is  estimated  to 
be  30  miles.  Although  the  standard  measure  of  variability  is 
the  standard  deviation,  all  error  estimates  made  in  the 
search  planning  process  are  In  terms  of  probable  error.  It 
will  be  demonstrated  that  this  error  Is  defined  by  the 
statistical  term,  "probable  error".  Survivors  are  believed 
to  be  In  rubber  rafts. 

The  proper   Coast  Guard  authorities  are   notified  at 


2100Z  on  22  June,  and  a  search  plane  has  been  dispatched.  It 
Is  due  to  arrive  on  the  scene  at  0100Z  23  June.  The 
aircraft's  navigation  error  (again  a  probable  error)  is 
estimated  to  be  10  miles. 

2.  Gathering  Data 

The  first  task  of  the  search  planner  is  to  estimate 
wind,  wind  current,  and  sea  current  at  the  position  of  the 
sinking.  Wind  Is  obtained  from  observations,  weather  maps, 
Navy  Fleet  weather  services,  or  the  U.S.  Weather  Bureau, 
whichever  is  available.  Wind  current  is  computed  from  the 
wind  by  means  of  a  graph  given  in  the  National  Search  and 
Rescue  Manual.  Average  Sea  Current  is  obtained  from  ocean 
current  charts  or  pilot  charts.  These  sources  for  average 
sea  current  are  compiled  from  observations  made  over  several 
years.  (Thus  this  data  may  not  be  accurate  for  the  time  and 
place  of  concern,  but  so  far  they  have  been  the  only  data 
available  to  the  search  planner.) 

3.  Drift  Plotting 
a.   Leeway 

Leeway  is  the  motion  of  the  survivor  craft  caused  by 
wind.  For  the  example  case  of  a  llferaft,  two  graphs  are 
provided  In  the  SAR  Manual  for  computing  leeway.  One  graph 
is  for  when  a  drogue  is  attached  to  the  raft  (a  drogue  is  a 
type  of  sea  anchor  designed  to  retard  liferaft  leeway)  and 
the  other  is  for  when  a  drogue  is  not  attached.  If  it  Is  not 
known   whether  a   drogue  is   attached,   both  conditions   are 


considered  separately. 

b.   Total  Drift 

Once  leeway,  wind  current,  and  average  sea  scurrent 
have  been  determined,  they  are  added  vectorial ly  to  produce 
total  liferaft  drift. 


T  c  r  ''1 L  C?fc»FT 


This  calculation  can  be  done  for  intervals  of  any 
length  -  one  hour,  three  hours,  twelve  hours  -  but  the 
interval  most  commonly  used  is  the  interval  between  wind 
observations.  That  is  if  wind  is  available  for  0900Z,  1200Z, 
1500Z,  etc.,  the  search  planner  would  probably  compute  total 
drift  in  three  hour  intervals,  from  0900Z  on  the  22nd  to 
0100Z  on  the  23rd.  Each  time  total  drift  is  calculated 
(another  probable  error),  12.51  of  that  distance  Is  retained 
as  "drift  error",  and  the  drift  errors  are  summed  to  obtain 
total  drift  error. 

Example:  In  this  manner,  the  position  of  the  example 
liferaft  is  calculated  to  be  U0-58.0N,  51-55. OW,  and  the 
drift  error  is  U.O  miles. 

c.   Minimax  Plotting 

If  it  is  unknown  whether  a  drogue  was  in  use, 
two  plots  are  made,  and  the  final  ("datum")  position  is 
defined  to  be  midway  between  the  two   positions  calculated. 
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The  drift  error  is  a  function  of  the  two  calculated  drift 
errors  and  the  distance  between  the  two  positions.  This 
process  of  plotting  two  positions  Is  called  MInimax  plotting 
i  n  the  SAR  Manual . 

d.   Error  of  Position 

Once  the  datum  has   been  computed/   it  is   then 
necessary  to  find  total  probable  error  of  position/  c,  where 

de  *  total  drift  error 

X   ■  error  in  initial  position  of  survivor  craft 

Y   ■  navigation  error  of  search  craft 


c  -  \l   de2+  XZ  -  Y* 

For  the  example/  c=  ^16+100+900  =^1016  *  31.9  miles. 

e.   Safety  Factor 

To  determine  the  size  of  the  area  to  be 
searched/  a  search  "radius"/  R  is  found  by  multiplying  c  by 
a  safety  factor/  depending  on  which  search  in  a  sequence  is 
being  planned.  For  searches  one  through  five/  the  respective 
safety  factors  are  1.1/  1.6/  2.0/  2.3/  and  2.5.  The  example 
of  the  sunken  fishing  boat  is  a  first  search/  so  R  ■  1.1  x 
1U.7  *  16.2  miles. 

The   area  to  be  searched/   then/  is   a  square  with 
sides  ■  2xR/  centered  on  datum: 
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f.   Sweep  Width 

The  standard  measure  of  the  effective  search 
width  of  a  search  craft  is  called  "sweep  width".  The  value 
of  the  sweep  width  is  obtained  from  a  table  in  the  Search 
and  Rescue  Manual  based  on  the  type  of  target,  search 
altitude,  and  meteorological  visibility.  The  sweep  width 
thus  obtained  is  corrected  for  the  effect  of  sea  state  by 
multiplying   by   a   "whitecap  correction   factor"  which   is 

tabulated  for  various  wind  speeds. 
C.   THE  COMPUTER  PROGRAM 

i.  objective 

The  computer   program  accomplishes  the  objective  of 
providing   an    estimate  of   the   survivor's    position   by 
automating  the  steps  of  the  manual  method  of  drift  plotting. 
2.   Source  of  Weather  Data 

In  choosing  a  source  of  weather  data  for  the  computer 
program,  three  criteria  were  considered:  a.  availability,  b. 
computer  access  time,  and  c.  accuracy.  The  final  selection 
for  data  source  was  the  Navy  Fleet  Numerical  Weather 
Central,  Monterey. 

a.   Availability 

FNWC  Weather   data  is   compiled  twice   daily,  at 
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OOOOZ  and  1200Z/  and  is  directly  available/  either  by  data 
link  or,  for  testing  purposes/  on  magnetic  tape.  For  use 
with  this  program/  the  data  is  stacked  chronologically  on  a 
master  tape. 

b.  Computer  Access  Time 

With  the  exception  of  the  mounting  of  the  master 
weather  tape,  all  weather  data  handling  is  eliminated. 
Access  speed/  then/  is  limited  only  by  the  speed  of  the  tape 
transport . 

c.  Accuracy 

Although  the  data  used  is  intended  for  wide  area 
coverage/  its  probable  error  is  less  than  the  probable  error 
assumed  when  computing  drift  error  (12.5%).  Although  there 
will  be  occasions  when  local  weather  reports  are  available 
and  may  be  more  accurate/  the  bulk  of  search  planning  is 
done  under  conditions  of  incomplete  weather  analysis  and 
forecasts/  and  when  gross  estimates  must  be  made.  Since  the 
weather  analyses  for  a  particular  time  period  and 
geographical  location/  once  assembled  for  this  program/  are 
fixed  parameters/  repeated  computations  of  the  same  search 
problem  will  always  produce  identical  answers.  This  permits 
a  rigorous  evaluation  of  the  data/  methods/  and  assumptions 
used  in  search  planning  by  a  systematic  testing  procedure. 
3.   Data  Description 

The   weather  data   are  produced   by   FNWC  for   3/969 
geographical  points  in  the   northern  hemisphere.   These  data 
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points  are  intersections  of  a  63  X  63  square  grid 
superimposed  on  a  polar  stereograph i c  map  projection.  Grid 
point  (32,32)  is  at  the  north  pole,  and  the  vertical  line 
( 1=32)  is  coincidental  with  the  80  W  and  100  E  meridians.  A 
sample  projection  is  shown  as  Figure  1. 

The  polar  stereograph i c  projection  provides  a 
plotting  reference  which  shows  an  entire  hemisphere  with  the 
least  possible  distortion  (Ref.  6),  and  the  scale  is  true  at 
60  degrees  North  Latitude.  Conversion  from  latitude  (L)  and 
longitude  (*)  to  grid  coordinates  is  accomplished  by  the 
following  formulas: 


(    =  32  ■+-   31,  20 


/-sin(J 


Cos  (a  t-io1) 


(l) 


--  yz  -  3i.  zc^ 


i  —  s  i  h 


L 


I  *~sih  L 


5tn  (X  1-  tOc) 


(2) 


For  the  example,  i  ■  38.008,  and  j  =  19.0143. 

The  data  is  stored  in  "fields"  of  3,969  words 
(numbers),  each  word  being  some  component  of  atmospheric  or 
oceanic  conditions  at  its  corresponding  (i,j)  grid  position. 
The  fields  required  for  direct  input  to  this  program  are: 


1U 


Figure  1.  Portion  of  Polar  Stereograph i c  Projection  Showing 
FNWC  (i,j)  Coordinates 
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U  Curr  U-  component  of  surface  current 

V  Curr  V-  component  of  surface  current 
U  Marine  U-  component  of  surface  wind 

V  Marine  V-  component  of  surface  wind 

where  the  U  direction  is  parallel  to  the  " i "  or 
horizontal  axis  and  the  V  direction  is  parallel  to  the  "j" 
or  vertical  axis,  in  the  FNWC  grid  system. 

The  Current  information  as  provided  by  FNWC  is  total 
surface  current,  including  wind  current. 

Fleet  Numerical  did  not  regularly  produce  U  Marine 
and  V  Marine  wind  fields,  so  these  fields  are  computed  by  a 
separate  computer  program  provided  by  the  Weather  Facility 
using  the  following  data  fields  as  input: 

D  1000         Height  of  the  1000-mb  pressure  line 
T  Air         Air  temperature  at  the  surface 
T  Sea         Sea  temperature  at  the  surface 

D.   CALCULATING  DRIFT  WITH  THE  COMPUTER 

l.  Leeway 

In  addition  to  the  liferaft  drift  tables  given  in 
the  Search  and  Rescue  manual,  there  are  other  ways  to 
predict  the  leeway  of  a  survivor  craft.  Most  of  these 
methods   involve   using   a   percent  of   the   wind   speed   to 
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determine  the  leeway  speed/  and  some  deflection  in  degrees 
to  determine  the  direction  of  leeway  off  the  wind  line.  For 
example/  a  boat  type  may  be  known  to  drift  in  a  direction  UO 
degrees  off  the  wind  line/  and  at  a  speed  equal  to  k%  of  the 
wind.  To  accomodate  these  methods  as  well  as  the  liferaft 
drift  table  from  the  SAR  Manual/  several  options  are 
available  to  the  user  of  this  computer  program. 

2 .  Options  for  Leeway 

a.  A  function  approximating  the  liferaft  drift 
tables  in  the  National  SAR  Manual  has  been  provided; 
drifting  speeds  for  wi  th  drogue  and  w?  thout  drogue  are  both 
aval labl e. 

b.  Leeway  speed  can  also  be  specified  as  a  percent 
of  wind/  in  whole  numbers/  from  2%    to  9%. 

c.  For  each  of  the  leeway  velocities  in  b.,  a 
deflection  off  the  wind  line  of  plus  or  minus  some  angle/  in 
tens  of  degrees/  from  0  to  90  degrees/  may  be  specified. 

d.  A  minimum-maximum  leeway  speed  option  with  the 
minimum  being  1%  of  the  wind/  and  the  maximum  being  some 
whole  number  percent  from  2%  to  9%  of  wind  speed/  is 
awa i 1 abl e. 

e.  No  leeway  at  all  (as  in  the  case  of  a  man  in 
the  water  or  some  other  almost  completely  submerged  object) 
may  be  specified. 

3.  Adjustment  of  Input  Data 
a.   Dates 


17 


All  dates  are  converted  from  month  and  day  to 
jullan  dates,  to  permit  reference  to  weather  data  and 
adjustment  of  the  program  "clock". 

b.  Latitude  and  Longitude 

All  initial  positions  are  converted  to  FNWC 
( i , j )  coordi  nates. 

c.  Leeway  Options 

Leeway  options  are  entered  by  means  of  a 
two-digit  code,  which  is  converted  to  speed  and  direction  of 
drift  within  the  program* 

d.  Weather 

Wind  and  current  are  provided  in  units  of 
cm. /sec.  and  are  scaled  in  the  program  to  nautical  miles  per 
hour.  They  are  further  scaled  to  decimal  fractions  of  the 
FNWC  grid  distance. 

k.     interpolation  pf  wither  Pat? 

a.   Interpolation  for  Position 

For  each  position  in  (i,j)  coordinates  there  are 
four  neighboring  grid  intersections  for  which  weather  is 
tabulated.   Example: 

GRID  POSITION:   i«38. 00787 

j=19.01U13 
GRID  POINTS:     .   (38,20)         .   (39,20) 

.   (38,19)         .   (39,19) 
All  weather  data  are  interpolated  linearly  among 
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these  four  neighboring  points.   For  example,  if  U  Current  at 

0000Z  at  grid  points  (i,j)  is: 

(38,20)  =  50.0     (39,20)  =  U9.0 
(38,19)  =  U2.0     (39,19)  =  UO.O 
then  for  i=38,  j=19.0H*13, 

U  Current  =  U2.0  +  .011*13  x  (50-U2)  =  1*2.112 
and  for  1=39,  a  similar  computation  results  in  U 

Current  =  1*0.12717. 

U  Current  for  the  position  is  then, 

U2.112  +  .00787  x  (-1.985)  =  U2.097  cm. /sec. 

b.  interpolation  for  Time 

For  each  hour  in  the  computation,  there  are  two 
observation  times  which  are  used  to  supply  weather  data: 
0000Z  and  1200Z,  one  of  which  will  be  prior  to  the  hour,  and 
one  which  will  be  after  the  hour.  A  linear  interpolation  for 
position  is  done  for  each  observation,  and  then  these  two 
values  are  interpolated  linearly  for  time.  For  example,  if 
the  U  Curr  for  1200Z  is  calculated  as  above  and  found  to  be 
37.25  cm. /sec,  the  interpolation  for  time  is  accomplished 
as  fol lows : 

U6.09  +  9/12  x  (-8.81*)  =  39.U6  cm. /sec. 

c.  The  Clock 

A   complete  set   of  interpolations   is  done   for 
each  position  and  for  each  weather  parameter  at  the  start  of 
the  drift  calculations,  and  each  position  is  advanced  for  an 
hour's   worth   of   drift.    This    interval   is   used   as   a 
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compromise  between  accuracy  requirements  and  computer  time. 
If  an  interval  much  smaller  than  one  hour  were  used,  a  large 
increase  in  computer  time  would  be  required.  A  problem 
which  would  take  20  seconds  at  one  hour  intervals  would  cost 
almost  20  minutes  of  computer  time  if  one  minute  intervals 
were  used. 

After  the  one  hour's  worth  of  drift  is  computed, 
and  the  positions  updated,  the  clock  is  advanced  one  hour, 
and  the  process  is  repeated  over  and  over  until  the  end  of 
the  problem. 

5.   Drift  Plotting 

Drift  calculations  are  applied  to  the  (i,j) 
coordinates  of  the  survivor  craft  plot.  That  Is,  the 
U-component  of  drift  Is  added  to  the  "1"  coordinate  of  the 
survivor  position  and  the  V-component  is  added  to  the  "j" 
coordinate. 

Example:  The  Interpolated  values  of  current  and 
wind  for  0900Z  on  22  June  are 

Parameter        cm. /sec.       miles/hour 
U  Curr         ♦  39.U60       ♦   .766 

V  Curr         ♦  36.128       ♦   .702 
U  Marine       +250.78        ♦  U.87 

V  Marine       -850.25        -16.55 

the  wind  speed  then  equals  v/4.8"?1*  I*-***  -17. 25 
nautical  miles/hour.  For  a  liferaft  without  drogue,  this 
wind  speed  produces  a  leeway  of   0.94  miles   per  hour  as 
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fol lows : 

Leeway  =  .226  xv/l7.25  =  0.94 

This  formula  is  an  approximation  to  the  liferaft 
drift  table  given  in  the  Search  and  Rescue  Manual. 

Since  arctan  ( 4. 87/-16. 55)  equals  163.6  degrees, 
this  angle  is  a  measure  of  the  wind  angle  (measured 
clockwise  from  the  vertical  (in  the  FNWC  system): 


If  leeway  deflection  is  specified  by  the  user/  it  is 
applied  to  this  wind  angle/  before  the  U  and  V  components  of 
leeway  are  calculated.  The  reduction  of  Leeway  and  wind 
angle  to  U  and  V  components  is  as  follows: 

U  Leeway  =  0.94  x  sin(163.6°)  =  0.266 

V  Leeway  =  0.94  x  cos(163.6°)  =  -0.902 

The  leeway  and   current  are  applied  to   the  position 
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as  fol lows : 

U  V 

Leeway        .266  -.902 

Current      . 766  . 702 

Total        1.032  -.200   knots 

These   values   are   then   divided   by   the   distance 
between  gridpoints,   at  this  latitude  ■  235.37  miles: 

.001*39         -.00085 
+  position   38.00787        19.01U13 
New  Position  38.01226        19.01328   for  1000Z  22  June. 
This  process  is  repeated  hour  by  hour  until  the  last 
datum   has   been   calculated/   drift   plotting   ceases,   and 
calculations  for  position  error  and  Search  Radius  begin. 
6.   Position  Error 

Just   as   in   the  manual   method,   the   initial 
position  error  of  the  survivor  craft  is  computed  as  follows: 
x  ■  distance  of  initial   (or  distress)  position  from 
last  known  position. 

A  «  error  of  last  known  position 

B  =  navigation  error  of   survivor  craft  as  a  percent 
of  distance  traveled 

The  initial  position  error  ,  X,  is 
X  »  A  ♦  B  x 

Since  Y   is  the  navigation   error   of   the  searcher, 
and  de  is   the  drift  error,  the  Total  Probable   Error,  c,  of 
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Position    is 

c    =    \7x*  ♦   Y*  +    dea 

7.  Safety    Factor 

c,  Total  Probable  Error  of  Position  is  multiplied  by 
the  safety  factor  to  find  the  Search  Radius  as  prescribed  in 
the  Search  and  Rescue  Manual. 

8.  Sweep  Width 

A   numerical    approximation    to   the   data    in 
Richardson's  Sweep  Width  tables  (Ref.  k)      may  be  seen  to  be: 
W   =  -3.3U377   +  1.25887   log  (1)   +  2.55UU2   log  (V)   + 
.01023  h 

where  1  is  boat  length 

V  is  visibility  in  miles,  and 
h  is  search  altitude  In  thousands  of  feet. 
In   addition,    corrections   to   Sweep  Width   for 
whitecaps  based  on  wind  speed  (v)  and  cloud  cover  (C)  can  be 
approximated  by: 

Fl  =  .886  +  .298X  -  .116  X2  (whitecaps,  X  =  v/10) 
F2  =  1.131  -  .U72C  +  .152C2  (cloud  cover) 

These  correction  factors  are  multiplied  by  W  to 
obtain  final  Sweep  Width.  This  Sweep  Width  function  is 
applicable  to  small  boats  only.  No  Sweep  Width  for  large 
ships  or  for  liferafts  is  available  in  the  program. 

E.   VERSIONS  OF  THE  PROGRAM 

The  Search  Planning  Program  is  written  in  FORTRAN  IV  for 
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use  on  a  remote  computer  terminal.  A  second  version  is 
written  for  use  on  the  operations  computer  at  Fleet 
Numerical  Weather  Facility.  No  remote  terminal  operation  is 
available  there,  so  these  two  versions  are  different  in 
several  ways: 

1.  Incut-Qutput 

a.  The  time  sharing  version  is  written  for  entry  of 
data  by  the  search  planner  onto  a  remote  console,  and  the 
results  are  printed  out  immediately  on  the  same  console. 
This  permits  the  planner  to  control  the  program  during 
execution,  and  to  make  intermediate  decisions  concerning  the 
problem.  An  example  of  the  output  of  the  time-sharing 
version  is  shown  on  the  next  two  pages.  Note  that  items 
typed  in  by  the  user  are  shown  in  lower  case,  and  the 
replies  typed  on  the  console  by  the  computer  are  in  upper 
case.  The  final  output  includes  sweep  width,  total  probable 
error  of  position,  search  radius,  and  calculated  wind. 

b.  The  batch-processing  version  at  FNWC  differs 
from  this  in  that  it  requires  cards  to  be  punched  for  input, 
and  the  output  Is  printed  in  message  form,  ready  for 
transmission  back  to  the  requestor. 

2.  Features  Omitted 

For  the  operational  version  of  the  program  at  Fleet 
Numerical  Weather  Facility  the  following  items  of  output 
shown  on  the  sample  time-sharing  printout  were  not  made 
avai lable: 
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$  sar 

EXECUTION  BEGINS. . . 

MONTEREY  SEARCH  PLANNING  PROGRAM 
ENTER  LEEWAY  OOOF,  1  RLANK,  THEN  DESCRIPTION 

>  11  1  iferaf t 

HOW  MANY  STARTING  POINTS  WILL  BF  ENTERED? 

>  1 

LIST  CHRONOLOGICALLY:  STARTING  TIM^S  AND  POSITION,  EX: 
0  70930Z  MAR, 3  8 2 2. 5 N, OB  ROB. 5 W 

>  220900z  jun, 4042. On, 05213.  Ow 

(+0       42.0     52       13.0    CHECK  :N    W 
CHECK.  TYPE  Y  OR  N 

>  y 

ENTER  NUMBER  OE  TIMES  DATUM  IS  TO  BE  COMPUTED  MOW: 

>  2 

ENTER  THESE  DTGS.   EX:020900  MAR, 031030  MAR 

>  230100  jun, 230700  jun 

SEARCH  1,  TRACK  1,  DATUM:  40  50.21  51  51.96 
SEARCH  2,  TRACK  1,  DATUM:  40  53.26  51  44.04 
ENTER  LAST  KNOWN  POSITION.  EX:  3522 . ON, 076 53 . OW 

>  MOO. On, 05400. Ow 
IS  LIFERAFT 

A  (l)BOAT,  (2)SHIP,  OR  ( 3 ) A  I RCRAFT?( ENTER  l,2,OR3) 

>  1 

ENTER  SEARCH  UNIT  FIX  ERROR,  Y.  EX:  5.0 

>  5.0 

INDICATED  DISTRESS  CRAFT  POSITION  ERROR, X: 

X=  15.00(1)  +   0.15(2)  TIMES  DIST.  FROM  LAST  KNOWN  POSITION. 

ANY  CORRECTIONS?  (TYPE  0,1,2,  OR  12) 

>  1 

ENTER  PARAMETER(l) : 

>  20.0 

IS  EACH  DATUM  TIME  A  DIFFERENT  SEARCH  ( 1, 2 , 3 , ETC . ) ? 

>  no 

ENTER  SEARCH  NUMBER  CORRESPONDING  TO  DATUM  TIME: 
**  1  0  SEARCH: 

>  1 

**  7  0  SEARCH: 
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>  o 

THF  FOLLOWING  ARF  THE  DATUM  POINTS  FOR  T"F 

1:   0  SFARfH: 
1*0-50. 2N      51-5?. OW     POS.  ERROR:   1*1.8,  RADIUS:   46,0  MILES; 

WIND  34.4KTS   FROM  267. 
SWFFP  WIDTH?   Y  OR  N 

>  y 

FNTER  SFARCH  ALTITUDE  (FX:800.) 

>  1000. 

ENTER  VISIBILITY  IN  MILES  (FX:10.0) 

>  20.0 

ENTER  BOAT  LENGTH  IN  FEET   (EX:36.0) 

>  10. 

FNTER  CLOUD  COVER  IN  FRACTION  (EX:  0.9) 

>  0.1 

SWEEP  WIDTH=       i*.6  Ml  LES. 

THE  FOLLOWING  ARE  THE  DATUM  POINTS  FOR  THE 

7:   0  SEARCH: 
40-53. 3N      51-UU. OW     POS.  ERROR:   1*1.9,  RADIUS:    0.0  MILES; 

WIND  34.4KTS   FROM  26  7. 
GOOD  LUCK  ON  YOUR  SEARCH!  LET  ME  KNOW  IF  YOU  FIND 
THELIFERAFT 


IHC002I  STOP       0 

R;  T=1.00/2.87  13.11.18 


Note:  Inputs  marked  with  ">' 
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a.  Printout  of  Intermediate  track  positions  (SEARCH 
1,  TRACK  1/  etc.)  Is  not  produced. 

b.  On  scene  Wind 

Since  the  wind   printout  as  shown  in   the  sample 
time-sharing  output  is  only  the  wind  velocity  and  direction 
at    the   beginning   of   the   search,   and   may   not    be 
representative  of  the  average  weather  conditions  throughout 
the  search,  this  value  is  not  included. 

c.  Sweep  Width 

Since  the  the  computation  of  sweep  width 
requires  data  which  may  not  be  immediately  available  to  the 
search  planner  at  the  time  of  his  request,  and  since  the 
formula  applies  only  to  small  boats,  this  feature  was  not 
included  in  the  Fleet  Numerical  version. 

d.  Leeway 

For  the  Fleet  Numerical  version,  leeway  is 
calculated  as  described  in  D.2.  The  time  sharing  version 
uses  a  square  root  function  for  all  leeway.  If  N  is  the 
percentage  drift  specified  (1%  to  9%),  the  leeway  function 
in  the  time  sharing  version  Is 

Ueway  =  i[5\J   WIND  (3) 

This  function  produces  leeway  which  Is  equal  to 
N|  of  the  wind  for  wind  speeds  of  16.0  knots.  For  other 
wind  velocities  the  function  follows  the  general  contour  of 
the  llferaft  leeway  (without  drogue)  table  in  the  Search  and 
Rescue  Manual . 
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III.  SEARCH  THEORY 

In  part  III  the  method  of  estimating  the  position  of  the 
search  object  and  the  rules  for  determining  the  size  of  the 
search  area  were  discussed.  This  is  only  part  of  the  search 
planning  process;  this  section  will  consider  how  resources 
are  assigned  to  the  search/  and  how  the  effectiveness  of  the 
search  is  measured. 

In  a  search  pattern  employing  parallel  sweeps,  the 
distance  between  these  sweeps  is  called  track  spacing,  S. 
The  Search  and  Rescue  Manual  recommends  that  for  urgent 
cases  the  track  spacing  used  be  equal  to  the  Sweep  Width,  W, 
and  for  less  urgent  cases  It  be  equal  to  twice  the  Sweep 
Width.  Once  the  value  of  track  spacing  has  been  decided 
upon,  the  amount  of  searching  required  (in  miles  flown)  L 
can  be  computed  from  S  and  the  area  A  to  be  searched: 

L  »  A/S 

The  measure  of  effectiveness  of   a  search  which  is  given 
in  the  Search  and   Rescue  Manual   is  called   Probability  of 
Detection,  and  it  is  presented   as  a  series  of  curves  plotted 
against  "coverage   factor"  or  W/S.   These  curves   are  called 
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"first  search",  "second  search",  etc.,  probabilities.  If  p 
is  the  probability  of  detection  in  the  first  search  curve, 
the  probabilities  2,..n..5   are  given  by 


i  -  (i  -  P)m 


These  curves  do  not  take  Into  account  the  size  of  the 
area  searched,  however,  so  that  a  search  can  be  conducted  in 
a  very  small  area  using  a  high  coverage  factor  which  will 
yield  a  very  high  probability  of  detection  according  to  the 
graph.  Common  sense  indicates  that  to  confine  the  search  to 
a  very  small  area  will  probably  not  result  in  the  detection 
of  the  target.  Nowhere  does  the  Search  and  Rescue  Manual 
define  the  probability  that  the  target  is  in  the  search 
area. 

The  absence  of  the  treatment  of  this  seemingly  critical 
factor  led  the  author  to  a  review  of  the  field  of  search 
theory  as  it  applies  to  survivor  search. 

The  basis  for  most  of  the  search  planning  criteria  in 
the  National  Search  and  Rescue  Manual  was  found  in  Koopman's 
(19UU)  work,  Search  and  Screening  (Ref.  2).  This  treatise 
developed  the  concept  of  sweep  width,  and  found  the 
probability  of  detection  for  various  situations. 

A.   DETECTION  LAWS 

1.   Instantaneous  Probability  of  Detection 
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If  Y  i  s  defined  so  that   the  probability  of  detecting 

a  target   in  a  short  period  of   time  dt   is  Ydt,   then  the 

quantity  "Y  is  called  the  instantaneous   probability  density 

of  detection.  When  the  searching  is   done  during  a   time  t 

under  constant  conditions,  the  probability  p(t)  of  detection 
is  given  by 


(t)=   /  - 


-yt 


f 

Search  and  Screening   developed   this   quantity   for 
various  situations. 

a.  Definite  Range  Law 

This  law  assumes  that  there  is  a  certain  range  R 
from  the  observer  beyond  which  It  is  impossible  to  detect 
the  target;  inside  this  range  detection  is  certain. 
Instantaneous  probability  of  detection  by  this  law  is 

y  (r)  -  oo      r  <  ft 

(5) 
t(0  =  0      r>R 

b.  Inverse  Cube  Law 

The  inverse  cube  law  is  the  result  of  making  the 
following  assumptions: 

(1)  The  observer  is  at  an  altitude  h  above  the 
ocean. 

(2)  The  observer  detects  the  target  by  seeing 
Its  wake. 

(3)  The  instantaneous   probability  of  detection 
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Is  proportional  to  the  solid  angle  subtended  at  the  point  of 
observation  by  the  wake. 

The  calculation  of   the  solid  angle  Is   shown  In 
figure  2. 


Figure  2.   Calculation  of  Solid  Angle 
Sol i  d  Angle  =  <x-  /3 

but  cA.  -  <i/5 

thus  t   ^o.h/s 


2.  \-7z 


sol i  d  angl e, 


Sol  i  d  Angle  ■   "-  ^  "  /  ' 

Since  y  was   assumed  to   be  proportional   to  the 

k  h 


where    k    included    all    factors    such    as 
meteorological  visibility,  etc. 
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Assuming  further  that  in  the  majority  of  cases  r 
is  much  larger  than  h,  then 

If  r"  ~ 

2.   Lateral  Range 

When  the  observer  and  the  target  are  on  their 
straight  courses  at  constant  speeds  for  a  long  time  before 
and  after  their  closest  point  of  approach  the  probability  of 
detection  may  be  shown  to  be  a  function  only  of  the  distance 
at  closest  point  of  approach.  This  distance  is  called  the 
1 ateral  range,  and  is  denoted  as  x.  This  probability  of 
detection/  p(x)/  plotted  against  lateral  range/  is  called 
the  lateral  ran, fie,  cyryg. 

Twice  the  area  under  the  lateral  range  curve  Is 
defined  as  effective  search  width  (or  Sweep  Width)/ 


T  Or> 


W-  */     f 


(x)d. 


(6) 


Koopman  shows  that: 

"The  effective  search  width  is  twice  the 
range  of  a  definite  range  law  of  detection 
which  Is  equivalent  to  the  given  law  of 
detection  in  the  sense  that  each  of  the  two 
laws  detects  the  same  number  of  uniformly 
distributed  targets  of  identical  velocity." 

Thus/  for  the  definite  range  law  of  detection  p(x)  » 

1  if  xl  R,  and  p(x)  »  0  if  x   >  R,    so  that  W  »  2R.  To  find  W 

for  the  inverse  cube  law/  first  p(x)  may  be  shown  to  be 
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p(x)=,    I    -    c 


2  kh 

GOV 


(7) 


where   w    is    searcher's    speed. 


From    (G)/ 


w  = 


I 


x  ^g 


2  M 


—       € 


)         ^\ 


y. 


«        2- 


zrrk  f\ 
cu 


(8) 
(9) 


3.   Random  Search 

To  find  the  probability  of  detecting  a  target  known 
to  be  in  an  area  A,  where  both  searcher  and  target  are 
moving  over  the  ocean  in  complicated  paths,  the  situation 
called  random  search  may  be  considered.  This  requires  the 
following  assumptions: 

a.  The  target's  position  is  uniformly  distributed 
in  A,  given  it  has  not  been  detected 

b.  The  observer's  path  is  random  in  A  in  the  sense 
that  it  can  be  thought  of  as  having  its  different  (not  too 
near)  portions  placed  independently  of  one  another  in  A. 

c.  W  is  small  with  respect  to  the  area  A. 

If  the  searcher's  path,  of  length  L,  can  be 
considered  to  be  composed  of  n  separate  and  equal  portions 
of  length  L/n,  then  the  probability  that  the  target  be 
detected  in  any  one  path  is  WL/nA.   The  probability  that  the 
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target  not  be  detected  at  all,  then  is 

.  n 


(10) 


and  for  large  n, 

-  W  l/A 

p  =  |  -  e  (id 

This  Is  the  formula  of  random  search. 
k.       Parallel  Sweeps 

The  most  common  method  of  non-random  searching  is 
one  which  uses  a  series  of  straight  parallel  sweeps  a 
distance  S  (or  track  spacing)  apart/  where  the  position  of 
the  target  is  fixed  with  relation  to  the  sweeps.  Let  the 
target  be  uniformly  distributed  In  an  area  of  size  A. 
Search  and  Screening  gives  the  following  result  for  P(L), 
the  probability  of  detecting  with  path  length  L=A/S. 

(1)  Definite  Range  Law 

Since  detection  will  surely   occur  if  the  target 

is  within  the   detection  range  W/2,  then   the  probability  of 

detection  Is  equal  to  the  coverage  factor  W/S. 

Thus, 
y>=  vVl_/A      l_1  A/xa/ 

p  =    I  L  >  A/W 

(2)  Inverse  Cube  Law 


(12) 


(13) 


3U 


(3)   Uniform  Random  Search 

From  (U), 
p=  |  —  e 

or  when  S^L/A, 


tV-/  A 


-  IA//5 


f  =- 


-e 


The   three  probabilities   using  parallel 
are    shown  in  Figure  3. 


(HO 

sweeps 


Pis) 


DEFINITE    RAN^E    l-AW 


RANDOM    5E.AR.CH 


INVE-RS£    cu&e   LAW 


W/5 

Figure  3.   Probability  of  Detection 

B.   DISTRIBUTION  OF  SEARCH  EFFORT 

l.  General 

Another  result  in  Search  and  Screening  involves  the 
definition  of  search  effort  density  as  a  function  of  (x,y) 
in  A. 

If  B  is  the  area  of  a   subreglon  in  A,  and  L   is  the 
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length  of  the  observer's  path  in  B,  then  the  expression  L  /B 

is  defined   to  be  the   observer's   track  density  in  B.    As  B 

becomes  very   small  (approaches   a  point  (x,y))  the   search 

densi ty  ft (x, y)  is  the  limit  of  WL  /B,  or  Wl(x,y),  where 

Jf  ((XjyJtxdy    =    L_  (15) 


and 


j(  (f  (x;  y-)dxrfy 


dyj  «    W/_  =£  (16) 


A 
Thus  if   (x,y)  is  applied   at  (x,y),  the  probability 

of  detecting  a  target  located  there  is 

^  -  /   e 

(17) 
Given   the  search   density  function  J^"(x,y)  and   the 
target's  distribution   p(x,y),  the  probability   of  detecting 

the  target  is  ^  (p(/  y) 

pO?)   -  f(p(*;Y)(l—e  )d*dy  (18) 

A 
The   fZ>  (x,y)   which   maximizes  P(  (f   )   is   shown   by 

Koopman  to  be 

r  (*j  y) 

(19) 


<p(x,y)=    log 


/ 
when  (x,y)  is  in  A/  and 

<p(*,  '/)    =     °  (20) 

when  (x,y)  is  not  in  A  , 

where  A  denotes  that  portion  of  A  in  which  the 
expression  (  <9  )  is  non-negative.  The  value  of  b  is 
determined  by  using  equation  (16). 
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2.   Normal  Target  Distribution 

If  the  target's  distribution  is  a  bivariate  circular 
normal   distribution   about  a   point(0/0)   in   a    cartesian 
coordinate  system, 


p(hy)  = 


r- 


2.7rcrz 


X       1       i. 

x    -  a    +■  y 


(21) 


then  the  search  effort  density  which  maximizes  P(  (p ) 


is 


<P<*,y)  =  -4f^" 


2  cr 


\f   <    CL 


(22) 


and 


4Ux,yJ  -  0 


r  ^o- 


(23) 


Thus   letting  a   be   the  radius   r   at  which   p(x,y) 
assumes  the  critical  value  b, 

/\^  -  TT^  ,  b  =  j^  e  ^*  (2%) 

where 


0.  - 


'4^  E 
7T 


i 

1 


(25) 


The  form  of  this  equation  is  an  inverted  parabola  of 
revolution  with  its  center  at  (0,0)  with  a  maximum  value  of 
ax/2c*. 
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IV.   THE  UNIFORM  COVERAGE  SQUARE  SEARCH 

Although  (  19  )  gives  the  "optimum"  search  effort 
di str i but ion/  this  optimum  distribution  is  difficult  to 
achieve  in  practice.  Flying  in  a  large  circular  area  with  a 
continuously  varying  effort  is  no  simple  navigational  task. 
Much  easier  is  the  common  practice  of  (1)  defining  a  square 
area  with  its  center  at  the  center  of  the  target's 
distribution  and  (2)  conducting  a  regular  search  of  uniform 
density  throughout  this  area.  If  the  target  is  not  found, 
then  (3)  the  area  is  expanded  and  the  search  repeated. 
Although  a  single  search  of  this  type  does  not  conform  to 
the  parabolic  optimum  search  effort  distribution,  a  series 
of  them,  each  larger  than  the  one  before,  may  be  a  very  good 
approximation. 

The  most  Important  advantage  of  this  plan  is  its  ease 
of  navigation,  allowing  the  pilot(s)  free  time  to  look  for 
the  target. 

A.   THE  PROBLEM 

The  problem  was  to  determine  the  optimum  size  of  a 
square  area  within  which  a  uniform  search  with  parallel 
sweeps  is  to  be  conducted. 
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B.   TERMS 

The  terms  and  definitions  in  the  Search  and  Rescue 
Manual  differ  from  those  In  Search  and  Screening.  In  the 
Search  and  Rescue  Manual/  probability  of  detection  is  used 
only  for  the  probability  of  detection  with  parallel  sweeps/ 
(given  a  uniform  target  distribution  within  the  search 
area).  In  order  to  avoid  confusion/  the  following  terms  are 
def i  ned: 


P  -  Probability  of  Success 


p  -  Location  Probability 


The  probability  that  on  a 
given  search  the  target  will 
be  found. 

The  probability  that  a  target 
is  within  a  specified  area. 


d   -   Probability  of  Detection    -    The   probability  of 

detecting  a  uniformly 
distributed  target  given 
paral 1  el  sweeps. 


W  -   Sweep  Width 


S  -  Track   Spacing 


L  -  Path  Length 


The  effective  search  width  of  the 
searcher;  the  area  under  the 
lateral  range  curve. 

The  mean  distance  between  adjacent 
parallel  sweeps  of  the 
searcher. 

The  total  number  of  miles  flown  by  the 
searcher  within  the  search 
area. 


c  -  Probable   Error  of  Position   -   Probable   error  of  the 

target's  position   from  datum; 
given   a   normal   distribution/ 
ca1.18  x  standard  deviation  of 
the  target's  distribution. 


R   -  Search   Radius 


The   search   radius   is  defined   as 
one-half  the  side   of  a  square 
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search  area. 


C.   THE  SOLUTION 

1.   The  Objective 

The  objective   is   to  maximize  the   Probability   of 
Success   P,  where   the  detection   function   for  a   uniformly 
distributed  target   known  to   be  in  the   area  is   defined  as 
d(W, L, R,x, y) .   P  is  expressed  as 


For  a  uniform  search,  the  probability  of  detection  is 
constant  and  equal  to  d(W,L,R)  inside  the  search  area  A,  and 
is  zero  outside/  (given  uniform  visibility  and  other 
factors),  and  so  d(W,L,R)  can  come  out  of  the  integral: 


p(WjL;c;e)»   cRWjM)      p(*,y)d*dy 


It  is  convenient  to  choose  x  and  y  as  distances  In  a 
cartesian   coordinate  system  with  the   target's  most   likely 
position  (datum)  at  (0/0). 

2.   Calculation  of  Individual  Terms 

Based  on  the   inverse  cube  law  for   parallel  sweeps, 
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Koopman  derives  the  following   expression  for  probability  of 
detection,  from  (i3): 

where  S  =  A/L  «  (2R  //L. 
This  is  equivalent  to  the  "First  Search"  probability  of 
detection  given  In  the  National   Search  and  Rescue  Manual. 

The  location  probability/  p,  is  given  by  the  formula 
for  the  circular  bivariate  normal  distribution: 


p 


(^ . 


I 


z<r 


ZTW 


i.y*  -* 


ax  a 


Where  R  defines   a  square  centered  at   (0,0)  with  sides 


■2R. 


Figure  U  is  a  graph  of  this  function,  along  with  a 
similar  curve  showing  the  cumulative  probability  for  a 
ci  rcl  e  with  radius  R.  Since  the  square  area  is  actually 
larger  (it  is  circumscribed  about  the  circle)  its 
probability  for  any  radius  R  Is  greater.  However,  If  these 
cumulative  probabilities  are  plotted  against  the  area  of 
the  search,  the  probability  for  a  square  will  be  slightly 
less  than  for  a  circle  (see  figure  5).  Although  there  is 
some  loss  of  efficiency  in  searching  a  square  instead  of  a 
circle,  this  loss  Is  very  minimal,  about  1%. 
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NOTE:  Although  the  standard  deviation  is  the  commonly 
accepted  measure  of  variability,  the  National  Search  and 
Rescue  Manual  uses  the  term,  "Total  Probable  Error  of 
Position".  Its  value  is  computed  using  estimates  for  drift 
error,  target's  position  error,  and  search  craft  navigation 
error.  These  errors  are  squared  and  added,  and  then  a 
square  root  is  taken  to  produce  Total  Probable  Error.  This 
term  Is  not  specifically  defined  in  that  publication,  so 
that  one  of  the  objectives  of  this  solution  was  to 
demonstrate  that  it  is  equivalent  to  the  statistically 
well-defined  term  'Probable  Error',  which  is  a  median  error: 
that  i S/  half  of  all  errors  have  absolute  values  greater 
than  that  error,  and  half  have  absolute  values  which  are 
smaller.  The  absence  of  this  definition  is  a  gap  in  the 
present  search  planning  process,  because  planners  are  called 
upon  to  make  estimates  of  errors  without  knowing  how  these 
errors  are  defined. 

3.  The  fisneral  Solution 

To  determine  the  optimum  size  for  a  square  search, 
given  a  certain  amount  of  searching  miles  available  (L), 
sweep  width  (W),  and  probability  distribution  (p(x,y)),  a 
straightforward  maximization  method  is  employed  In  the 
following  sequence: 

a.  Determine  total  effort,  E. 
E=WxL 
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b.  Choose  some  arbitrary  initial  search  radius/  x. 

c.  Determine  the  search  area,  A,  for  this  x: 
A  =  4  x% 

d.  Determine  effort  density  for  a  uniform  search: 

e.  Determine  probability  of  detection  d  from  ( i3  ) . 

f.  Determine  the   probability  p  from  (2.1)   with  R  = 


x. 


g.   Determine  the  probability  of  success, 
P(x)  =  p(x)d(x). 

h.  Increase  x  by  Ax  and  repeat  calculations.  If 
P(x)  increases,  continue  adding  Ax  until  a  local  maximum  is 
encountered.  If  P(x)  decreases  initially,  reduce  x  by  2Ax 
and  continue  subtract  i ng  Ax  until  a  local  maximum  is 
reached. 

Since  a  cumulative  probability  distribution  is 
continuous  and  has  a  slope  throughout  its  range  of  2.  0,  and 
since  the  probability  of  detection  function  d(x)  is 
continuous  with  a  slope  of  <  0,  the  function  P(x)  =  p(x)d(x) 
can  be  shown  to  be  unimodal,  and  so  a  local  maximum  in  P(x) 
will  be  an  absolute  maximum. 
4.   AdpI icat  ion 

The  first  use  of  this  method  was  to  check  the 
National  Search  and  Rescue  Manual's  recommended  search 
planning  guidelines  to  see  if  they  produce  square  searches 
of  optimum  sizes.   These  guidelines  are  given  in  the  form  of 
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'Safety  Factors1,  or  search  radii,  for  repeated  searches,   as 
multiples  of  Total  Probable  Error  of  Position,  e.g., 

Search  Safety  Factor 

1  1.1 

2  1.6 

3  2.0 
k  2.3 
5  2.5 

As  far  as  effort  density  is  concerned,  there  are  two 
recommendations  given  in  the  SAR  Manual.  One  is  for  a  track 
spacing  equal  to  sweep  width,  for  urgent  cases,  and  the 
other  for  a  track  spacing  equal  to  twice  the  sweep  width, 
for  less  urgent  searches.  Both  of  these  recommendations 
were  tested. 

The  approach  used  was: 

a.  Set  some  arbitrary  value  for  W 

b.  Set  S=W  so  that  WL/A.  =1. 

c.  Find   R   from  "safety   factor"   table   for   the 
appropriate  search  (starting  at  1). 

d.  Find  path  length,  L: 

-   4  ft* 
U  '    5 

e.  Beginning  at  x  ■  R,  and   with  W  and  L  above,  use 
the  general   solution  process  described   in  paragraph   3.  to 


U6 


find  x',  the  optimum  search  radius  for  that  path  length. 

f.  Assuming  that  the  opt  imal  search  is  executed 
unsuccessfully,  reevaluate  the  target's  distribution  p(x) 
(see  Para.  5.)  and  repeat  the  process  from  c.  above.  After 
five  successive  searches  have  been  computed,  stop  the 
computation. 

g.  Compare  the  resulting  five  optimum  sizes  with 
the  sizes  recommended  in  the  Safety  Factor  table  from  the 
Search  and  Rescue  Manual. 

5.   Baveslan  Analysis 

Bayes'  Law  of  Probability  says  that  an  estimate  of 
the  probability  of  a  hypothesis  can  and  should  be  revised  on 
the  basis  of  the  outcome  of  a  related  event.  In  the  case  of 
search  theory,  the  hypothesis  is  that  the  target  is  inside  a 
given  area.  A  search  is  conducted,  and  the  event  is  either 
the  detection  or  the  non-detection  of  the  target.  This 
event  should  revise  the  initial  estimate  of  the  probability 
that  the  target  was  in  the  area.  Of  course,  If  the  event  is 
detection,  the  new  probability  that  the  target  is  In  the 
area  Is  1.  If  the  event  is  non-detection,  then  Bayes' 
formula  must  be  used  to  calculate  the  new  probability. 

1.   Exampl e 

The  probability  that  a  target  is  in  a  certain  search 
area  has  been  determined  to  be  50%,  and  a  search  is 
conducted  in  this  area  with  a  uniform  probability  of 
detection  of  50$,  The  resulting  probability  of  success  is  .5 
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x  .5  =  .25  or  25%. 

For  the  use   of  Bayes'  formula,  the   following  terms 

are   defined: 

A    The   hypothesis  that  the   target  is   in  the 

area. 

A    The   hypothesis  that  the   target  is   not  in 

the  area. 

D     The  event  of  detection. 

D    The  event  of  non-detection. 


Then 


p(AlD)  = 


(A)  •  p  (pU) 


p  W  ■  plDWX)  *p(A)  •  p(DiA) 


When  the  target  is  detected: 


p(/»lD)  ^ 


0.5"  x     0  b" 


0.5"  *  OS      i-O* S 


0.0 


O 


as  expected.   If  the  target  is  not  detected,  the  new 
probability  would  be 

p(a)  •  p(5\a) 


f UId)  = 


K*)  '  f(D)b)     *  p(£  ;«  p(5|A) 


0<5"  /  0-5  ■»■  O.s-  x  /,0 


3 
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Whereas  the  target  had  a  50%  a_  priori  probability  of 
being  in  the  search  area,  its  new,  or  a.  poster  i  or  i 
probability  of  being  in  the  area  is  33.3%. 

For  the  solution  in  paragraph  U.,  the  entire 
probability  distribution  of  the  search  object  was 
reevaluated  between  successive  searches.  This  was  done  for 
p(x)  inside  the  search  area  by 


/(x)  ,     p(x)O-d) 


d    -      pU)d ) 


(26) 


and  for  p(x)  outside  the  search  area  by 

(>    -  P<-P-U    )  (27) 

where   R  is  the  radius  of  the  search  conducted. 


In  order  to  illustrate  the  general  form  of  the 
probability  distributions  between  searches,  a  function  was 
used  which  is  similar  to  a  density  function.  This  function, 
f(x)  was  defined  to  be 

dO(xJ) 

(28) 

where  x  is  the  measure  of  radius,  or  again,  one-half 
of  one   side  of  a   square,  and   8xdx  is  the   differential  of 


f^=  (^7—  )[ir) 


1(9 


area  ( F  i  gure  6  ) . 


-  8xdx 


Figure  6.  Derivation  of  Density  Function  for  Square 
Area 

Curve  1  in  Figure  7  is  the  density  function  before 
any  searching  is  done/  as  derived  by  (28)  from  (21).  Curves 
2  through  5  are  these  functions  after  successive  searches. 

Figure  8  is  a  similar  set  of  curves  showing  the 
cumulative  probability  distributions. 

6.  Results 

This  method  (from  k.)  gives  the  following  optimum 
searches: 
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Figure    7. 
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Probability  Density  with  Application  of  Bayes' 
Between  Successive  Optimum  Square  Searches. 
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U  2.U3  2.3 

5  2.71+  2.5 

Since  these  results  were  very  close  to  the  'Safety 
Factors'  listed  above/  It  was  concluded  that  (1)  the  safety 
factors  were  originally  good  cal cul at  ions,  and  (2)  they  were 
calculated  using  the  mathematical  definition  of  probable 
error. 

It  should  be  noted,  however,  that  these  safety 
factors  are  close  to  optimum  search  radii  only  for  the 
condition  where  track  spacing  equals  sweep  width.  The 
second  recommendation,  track  spacing  equals  twice  sweep 
width,  does  not  produce  anywhere  near  optimum  search  areas. 
When  a  track  spacing  S'  equal  to  twice  the  sweep  width  is 
used,  the  new  path  length  L'  is  given  by 

li!   _  iJL*    _L 

For  this  new  path  length,  optimum  search  sizes  were 
found  to  be: 

Search  Optimum  Radius 

1  0.85 

2  1.30 

3  1.66 
U  1.95 
5  2.20 
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When  there  Is  less  search  effort  available  the  size  of 
the  search  area  as  well  as  the  effort  density  should  be 
reduced.  The  SAR  Manual  gives  no  guidance  as  to  how  much  the 
area  should  be  reduced  with  a  given  amount  of  available 
effort.  Conversely/  there  Is  no  Indication  as  to  how  excess 
search  resources  should  be  used/  If  such  a  condition  should 
exist. 

What  is  needed  is  some  way  for  the  search  planner  to 
determine  the  optimum  search  configuration  directly  from  the 
amount  of  search  effort  available. 

D.   THE  GENERAL  CASE 

Assume  that  a  series  of  optimum  square  searches/  each 
larger  than  the  preceding  one/  is  a  good  approximation  of  an 
optimum  search  effort  distribution.  Consider  all  the 
searches  of  the  series  combined  Into  one  search  effort/  so 
that  the  search  effort  density  <Pir)  assumes  a  step  function 
which  is  close  to  the  ideal  effort  denstty. 


t 
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The  radius  r  of  the  combined  search  Is  equal  to  the 
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radius  of   the  last  (and   the  largest)  individual   search  in 

the  sequence. 

If  the  search  had  been  a  circular  one,  and  if  the  effort 

density  had   conformed  to   the  optimum,   the  expression   for 

this  radius  would  be 

1 
4cr2   £ 

a  - 


rr 


But  since  the  search  is  composed  of  a  small  number  of 
square  searches,  this  same  formula  is  not  expected  to  apply. 
There  should,  however,  be  some  way  to  approximate  a  function 
which  will  define  the  optimum  radius  based  on  available 
effort  and  the  target's  distribution.  The  simplest  way  to 
do  this  is  to  compute  total  search  effort  for  a  few  optimum 
square  searches,  plot  the  radius  of  each  of  the  searches 
against  the  cumulative  search  effort,  and  attempt  to  smooth 
a  curve  through  them. 

This  was  done  for  the  ten  searches  (five  at  W  =  S  and 
five  at  W  =  2S)  calculated  above,  and  the  resulting  curve  is 
shown  in  Figure  9.  Just  as  radius  is  in  terms  of  probable 
error,  so  are  the  factors  W  and  L  of  total  effort 
'normalized'  by  dividing  each  by  probable  error.  Thus 
'normalized'  effort  is  a  multiple  of  the  square  of  probable 
error . 

It   should   be   cautioned   that  the   curve   is   only   an 
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Figure  9.  Optimum  Size  of  Square  Search  based  on 
Cumulative  Normalized  Effort  Factor. 


0.5  i.o  I    5 

Search      Radius, 


2.0 
R  /   c 
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approximat i on/  and  will  not  be  valid  if  the  effort  available 
for  a  single  search  is  very  large.  Thus  if  a  search  is 
computed  which  shows  a  coverage  factor  (W/S)  greater  than 
1.8,  then  the  results  of  the  curve  will  not  give  an  optimum 
search,  and  some  other  method  will  be  required.  Again,  the 
effort  plotted  is  cumulative,  so  that  all  the  effort  exerted 
in  earlier  searches  must  be  added  to  the  effort  available 
for  the  present  search  before  entering  the  table. 

Exampl e:  One  150-knot  aircraft  is  available  for  six 
hours  of  searching;  the  total  probable  error  of  position  has 
been  calculated  to  be  31.9  miles,  and  W  is  estimated  at  10 
mi  1 es. 

Total  effort  then  Is  found  as  follows: 

E=  150kts  X  6hrs  X  lOmiles(W) 

.   1     GO 
-    9000  square  mi les 

In  terms  of  probable  error,  (c=31.9) 

ECnormal Ized)  =  9000/(31.9  X  31.9)  =  8.9 


From  the  graph,  r  =  1.4, 


R  =  l.i*  X  c  =  i*i*. 6  miles,  or  89.2  miles  on  a  side,  and 


S  =*  k    X  R  X  R  /  L   =  8.75  mi  les. 
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If  the  target  is  not  found  on  the  first  search,  and  a 
second  is  to  be  conducted,  the  amount  of  available  effort  is 
calculated  for  the  second  search.  The  normalized  effort  for 
this  search  must  be  combined  with  all  previously  expended 
effort  before  entering  the  graph.  Thus  if  E( norma  1 i zed )  for 
the  second  search  is  calculated  to  be  10.0,  the  graph  should 
be  entered   with  a  figure   of  10.0+8.9   or  18.9  to   find  the 

radius  of   the  second   search.    Each   search   area  will   be 
larger  than  the  preceding  one. 

This  decision  curve  should  prove  to  be  a  more  flexible 
tool  for  search  pi anni ng.  than  the  Safety  Factor  table,  since 
it  provides  a  continuous  range  of  search  radii  for  different 
amounts  of  effort  available. 
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V.   THE  SECTOR  SEARCH 


There  is  one  search  pattern  in  general  use  which  does 
effect  a  greater  concentration  of  effort  at  the  center  of 
the  search  area,  and  that  is  the  sector  search.  This  search 
pattern  is  executed  by  a  series  of  sweeps  crossing  over  the 
center  of  the  search  area  at  various  angles.  These  sweeps 
are  connected  at  their  endpoints  by  chords  flown  near  the 
perimeter  of  the  area.  This  pattern  covers  a  circular  area 
and,  perfectly  navigated,  produces  a  track  as  shown  in 
Figure  10. 

A.   SEARCH  EFFORT  DISTRIBUTION 

The   total  effort   expended  in   a  sector   search  can   be 
expressed  by  the  integral 

f  =  (     <p(r)      2*rr  dr 

where  R  is  search  radius 

r   is  the   distance   from   the  center   of   the 
search 
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Figure  10.  The  Sector  Search 
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(p  (r)   is  the   search  effort   density  at   some 
distance  r. 

This  density  function  @(r)  is  simply  W  times  the 
number  of  paths  which  intersect  a  circle  of  radius  r, 
multiplied  by  dr  to  give  total  path  length  between  r  and 
r+dr,  and  then  divided  by  the  area  between  r  and  r+dr.  Thus 


<jj(r)  = 


Then, 


2  IT  r    dr  2 IX  r  (19) 


C*£> 


Wn 


2nrdr      =      VN/nR. 


(30) 


and, 

L  =  nR 
as  expected. 
Here  n  is  the  number  of   'spokes'  produced,  or  twice  the 
number  of   complete  sweeps  across  the    circle.   This  density 
function  is  illustrated  in   Figure  11. 

B.   PATTERN  ERROR 

There  will   always  be  some   error  in  flying   any  search 
pattern,  which   should  be  taken  into   account   for  the  sector 
search.   The   navigation  error  of   the  searcher   has  already 
been  taken   Into  account  when   computing  the   Total  Probable 
Error   of  Position.    The   pattern   error,  however,   can   be 
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Figure  11.  Search  Effort  Density  In  a  Perfectly 
Navigated  Sector  Search. 
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different  from   this  navl gat  Ion   error.   Whereas   navigation 

error  considers  only  the  searcher's  ability  to  find  the 
datum  point  In  order  to  define  the  search  area/  pattern 
error,  as  It  Is  used  here,  describes  his  flying  error  within 
the  area.  If,  once  the  searcher  finds  what  he  thinks  to  be 
datum,  he  drops  a  smoke  float  or  some  other  marking  device, 
then  his  pattern  error  will  be  much  smaller  than  his 
navigation  error.  Conversely,  if  when  arriving  in  the 
search  area  he  stops  refering  to  his  navigation  aids  and 
relies  thenceforth  on  dead  reckoning  only,  his  pattern  error 
will  be  greater  than  his  navigation  error. 

Although  by  its  nature,  a  flying  error  is  bivarlate 
(like  the  error  of  target  position)  a  simplification  was 
made.  Since  the  pilot  Is  flying  long  sweeps  which  have 
their  endpoints  relatively  far  away  from  the  area  of  concern 
(the  center  of  the  area),  only  his  lateral  position  error 
(which  is  univariate)  was  considered. 

It  was  assumed  for  this  analysis  that  the  lateral 
component  of  pattern  error  is  nearly  normal,  and  that  It  can 
somehow  be  calculated  or  estimated. 

The  expression  for   the  total  effort  expended   In  flying 

the  spokes  of  a  sector  search  Is 

*  IE 
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E=   'zWn    7=—   e       r^rcl*  (31) 
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The  derivative  of  this  expression   with  respect  to  r  was 
used   as  the   search  effort   density  function   for  a   sector 


search  with  pattern  error. 
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For  the  analysis  this  basic  density  function  was 
computed  with  sweep  width,  standard  deviation  of  pattern 
error,  and  number  of  spokes  all  equal  to  1.  To  find  the 
density  function  for  any  particular  set  of  parameters,  the 
density  is  multiplied  by  W  and  n,  and  the  horizontal  scale 
is  multiplied  by  the  standard  deviation  of  the  pattern 
error. 

C.   OPTIMUM  SECTOR  SEARCH 

The  plot  of  this  density  function  is  found  in  Figure 
12.  An  inspection  of  this  plot  will  show  that  it  is  close  to 
a  quadratic  function  out  to  r  =  3  standard  deviations. 
Since  Koopman's  ideal  search  is  a  quadratic,  some  guidance 
was  provided  as  to  the  optimum  radius  for  the  sector  search, 
based  on  the  pattern  error.  It  was  then  necessary  to  adjust 
the  number  of  spokes  to  arrive  at  the  Ideal  effort 
expenditure  to  match  this  radius.  Koopman's  ideal  search 
effort  density  is  a  function  of  search  radius  and  target's 
pos  i  t  ion  error : 
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Figure  12.  Search  Effort  Density  in  a  Sector  Search 
with  Pattern  Error. 
65 


(f(',y) 


r  >  ^ 


O 
and  the  radius   is  a   function  of  target  position  error 


and  available  effort: 


a.  - 


-4^ 


IT 


1_ 
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Solving  for  effort/  E: 

a4  TT 


e  * 


4   cr 


W  L- 


However/  in  the  case  of  the  sector  search/  L,  the 
number  of  miles  flown,  is  simply  the  number  of  spokes,  r\, 
times  the  radius.  The  chords  flown  at  the  perimeter  are 
disregarded  in  this  analysis: 

L  =  a  X  n 


Substituting/  and  solving  for  n: 


a?  V 


4er*  W 

Ideally/  then/   since  the  optimum   a  is  known   to  be 


66 


close  to  three  standard  deviations: 


m 


2 1    <?f      TT 

4  cr2  W 


(33) 


or,  when  all  errors  are  expressed  as  probable  errors, 
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.67  O-k 


C 


\8  & 


cc\  ■= 


W    C1 


(3l») 


67 


Table  1  is  a  tabulation  of  this  function  for  various 
conditions  of  total  probable  error  of  position  c,  and 
pattern  error.  The  n  derived  from  the  table  must  be  divided 
by  sweep  width  W  to  obtain  number  of  spokes. 

From  the  example  of  the  sunken  fishing  boat/  c  ■  31.9 
and  W  =  10.0  miles.  If  the  searcher's  pattern  error  were  13 
miles,  then  search  radius  R  would  equal  58.5  miles. 

Interpolating  in  the  table,   n  =  100,  then 
Number  of  spokes  =  100/10  =  10 

and  the  angle  between  adjacent  search  legs  ■  360/10  =  36 
degrees. 

Total  search  miles  flown  would  be 
10  X  58.5  =  585 

plus  the  amount  of  flying  at   the  perimeter  (in  order  to 
fly   the  pattern,   the   aircraft   must  traverse   a   distance 
somewhat  less  than  one-half  the  circumference  of  the  area): 
585  +  58.5  X  TT  »  585  +  18U  =  770  miles. 

This  table  for  sector  search  will  not  be  useful  until  an 
accurate  means  of  predicting  pattern  error  Is  found.  The 
function  given  In  (34  )  Is  extremely  sensitive  to  pattern 
error,  and  thus  accuracy  for  this  parameter  Is  required  to 
use  the  table. 

Even  If  It  Is  not  actually  used  to  determine  the 
parameters   of  a   sector   search,  the   table   does  serve   to 
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illustrate  that  accurate  navigation  does  not  necessarily 
produce  the  optimum  sector  search/  and  that  only  under  the 
unusual  conditions  of  an  extremely  accurate  estimate  of 
target  position  is  it  desireable  to  reduce  pattern  error. 
On  the  contrary/  in  most  instances/  effectiveness  may  be 
enhanced  by  increas  i  ng  pattern  error. 
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VI .   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   THE  SEARCH  PLANNING  PROGRAM 

It  Is  recommended  that  the  Coast  Guard  begin  immediately 
to  make  full  use  of  the  Search  Planning  Program  now 
operational  and  undergoing  testing  at  the  Fleet  Numerical 
Weather  Central ,  Monterey.  A  guide  for  users  of  this  program 
is  included  as  Appendix  A. 

1.   Precautions 

The  following  cautions  are   appropriate  for  users  of 
the  information  provided  by  this  program. 

a.  Region  of  Validity. 

Because  the  weather  data  points  used  in  this 
program  are  approximately  200  miles  apart,  drift 
calculations  provided  for  areas  within  100  miles  of  shore 
should  be  double-checked  by  hand  before  using.  If  hand 
calculations  produce  a  larger  value  for  total  drift,  then 
the  hand-calculated  answers  should  be  used  in  lieu  of  the 
computer  solution. 

b.  Time  of  Datum 

Weather  predictions  in  the  format  required  for 
this  program  are  not  produced  by  the  Weather  Facility  for 
more  than   U8  hours  In   advance.   Accordingly,   users  should 
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avoid  requesting  plots  for  time  periods  beyond  36  hours  from 
the  time  the  request  is  submitted.  If  a  plot  is  requested 
which  is  later  in  time  than  the  last  computed  weather,  the 
information  will  be  provided  for  that  plot  and  a  caution 
statement  will  appear  in  the  message  reply. 

c.   Numbers  of  Initial  and  Datum  Positions 

No  more  than  nine  initial  positions  may  be 
entered,  and  no  more  than  nine  datum  times  will  be 
calculated,  In  any  one  run  of  the  program. 

The  user  Is  encouraged  to  request  whatever 
number  of  datum  points  may  be  useful  to  him.  Thus  if  a 
search  Is  planned  for  0800Z,  the  search  planner  may  wish  to 
request  additional  data  for  target  position  at  1000Z,  1200Z, 
UOOZ,  and  1600Z. 

2.  Leeway 

It  is  suggested  that  further  experimentation  be  done 
in  the  area  of  survivor  craft  leeway  to  determine  the 
general  shape  of  the  leeway  function.  It  is  expected  that 
this  shape  will  be  close  to  that  of  the  liferaft  drift 
tables,  and  to  the  formula  of  Witting  given  in  Ref.3. 


V  -   k  \j   Wind  (35) 

and  which  Is  used  In  the  time-sharing  version  of  the 
program  for  all  leeway.  This  same  function  is  used  in  both 
versions  to  approximate  liferaft  leeway. 
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When  better  data  is  available/  the  program  at  Fleet 
Numerical  can  be  changed  to  accomodate  the  curved  leeway 
function. 

3.  Testing 

It  is  suggested  that  a  rigorous  program  of  testing 
and  evaluation  of  the  results  of  the  Search  Planning  Program 
be  initiated.  If  the  program  proves  to  be  more  accurate 
than  other  methods/  Its  contribution  to  the  search  planning 
process/  and  to  oceanographic  research  In  general/  will  be 
substant 1 al . 

\*      Operational  Use  in  the  Coast  Guard 

Although  the  program  will  be  available  to  the  Coast 
Guard   via  the   facilities  of   the   Fleet  Numerical   Weather 

Central/  there  are  further   possible  advantages   of  such   a 
program  which  are  not  available  through  these  facilities. 

a.  Optimum  Search  Areas 

The  computer  method  of  solving  for  optimum 
searches  shown  in  part  IV  of  this  thesis  could  be  added  to 
the  Search  Planning  Program. 

b.  Other  Applications 

Computation  of  on-scene  aircraft  endurance/ 
time  enroute  search  areas,  division  of  search  area  among 
different  search  units,  maximum  rescue  coverage  intercepts/ 
fixes  from  radio  direction  finder  bearings,  and  many  other 
calculations  familiar  to  search  planners  could  all  be 
available  instantly  from  a  remote  computer  terminal. 
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The  Coast  Guard  Amver  computer  center  at  Governor's 
Island  New  York  regularly  receives  weather  data  by  wire  from 
Fleet  Numerical  Weather  Central  which  is  used  to  produce 
weather  maps.  The  fields  required  for  the  Search  Planning 
Program  could  also  be  transmitted  to  a  central  Coast  Guard 
computer  facility.  Remote  terminals  connected  to  this 
facility  and  installed  in  Rescue  Coordination  Centers 
(RCC's)  could  provide  the  output  from  this  program,  Amver 
surface  picture  plots,  and  other  useful  data  directly  to 
search  planners. 

3.   DETERMINATION  OF  SEARCH  AREAS 

It  is  recommended  that  the  findings  of  parts  IV  and  V  of 
this  thesis  be  evaluated  by  the  Coast  Guard  and  applicable 
portions  be  included  in  the  National  Search  and  Rescue 
Manual . 

C.   FUTURE  RESEARCH 

There  are  many  unexplored  areas  in  the  field  of  Search 
Planning  in  which  further  research  may  yield  valuable 
resul ts. 

1.   The  Trackl ine  Search 

Although  the  type  of  point-datum  search  presented  in 
parts  IV  and  V  are  computed  in  two  dimensions,  there  are 
times  when  the  target's  distribution  is  centered  on  a  line 
of   position   (or   trackline),    which   presents   a   problem 
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essentially  one  dimensional  in  nature. 
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If  a  boat  departs  point  A  intending  to  follow  a 
certain  track  AB  to  point  B,  and  does  not  arrive  at  point  B, 
and  if  no  other  information  is  available  as  to  his  probable 
position/  a  one-dimensional  or  trackline  search  problem  is 
presented. 

a.  Partial  Solution 

Assuming   that   it   is    decided   to   begin   the 
searching   efforts  by   executing  just   one   sweep  along   the 
length   of  the   trackline/  what  Is   the  probability   that  the 
target  will  be   located?   The  following  analysis   is  for  the 
inverse  cube  law  of  detection. 

Given   the  standard   deviation  of  the   target's 
error  of  position  CT ,    and  given  sweep  width/  VI,    where 


W-  */" 


zir  m 


(36) 


Solving   for  m,   which   includes  visibility   and 
other  factors/ 


m  = 


8  Tf 


(37) 
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The    probability   of    detecting   a    uniformly 
distributed  target  is 
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Introducing  a  new  variable,  /&  =  ~yj~ 
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The  probability   that  the   target  is   located  at 
some  distance  x  from  the  trackllne  may  be  assumed  to  be 
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Thus,  the   probability  of  detecting   this  target 
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which  has  the  solution 
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b.   General  Solution 

The  solution  for  probability  P  of  detecting  a 
target  on  a  trackline  search  for  more  than  one  sweep 
involves  a  more  complicated  expression,  one  which  does  not 
lend  itself  to  such  a  simple  solution.  Since  a  trackline 
search  is  executed  one  sweep  at  a  time,  a  general  solution 
for  the  optimum  distance  of  each  sweep  from  the  trackline  is 
Indicated,  and  should  be  the  subject  of  further  study. 
2.   Optimum  Search  Speed 

The  Search  and  Rescue  Manual  prescribes  a  "moderate" 
speed  for  searching.  Koopman's  formula  for  sweep  width,  from 
( 3   )  gives 
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This  indicates  that  sweep  width  is  inversely 
proportional  to  the  square  root  of  the  relative  speed  of  the 
target,  w.  When  search  speed  is  very  large  In  comparison  to 
target  speed,  w  takes  on  the  value  of  the  searcher's  speed. 
As  w  increases,  however,  the  amount  of  search  path  length  L 
that  can  be  covered  in  a  certain  period  of  time  increases 
proportionally: 
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It  follows  that  the  amount  of  search  effort/  W  x  L/ 
that  can  be  effected  by  a  certain  number  of  aircraft  in  a 
fixed  period  of  time  increases  as  the  square  root  of  search 
speed : 
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It  may  be  argued/  then/  that  when  time  is  a  factor/ 
such  as  where  persons  In  distress  may  not  survive  for  long 
periods  of  time/  searching  should  be  accomplished  at  the 
highest  speed  practicable.  In  addition/  since  the  coverage 
of  a  designated  search  area  will  take  less  time  at  a  higher 
speed/  lookouts  may  encounter  less  fatigue/  and  the 
effectiveness  of  the  search  may  be  enhanced. 

D.   CONCLUSION 

In  conclusion/  there  Is  a  need  for  more  scientific 
research  in  the  practical  aspects  of  search  planning,  and  a 
corresponding  need  for  the  education  of  search  planners  In 
the  fields  of  probability  and  search  theory. 
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APPENDIX  A 

USER'S  GUIDE  TO 
SEARCH  PLANNING  PROGRAM 


1.  Pescriptlon 

The  MONTEREY  SEARCH  PLANNING  PROGRAM  Is  a  computer 
program  which  solves  open  ocean  drift  plotting  problems/  and 
provides  datum,  minimax  datum.  Total  Probable  Error  of 
Position,  and  Search  Radius  for  both  simple  and  complex 
searches.  The  use  of  this  program  is  available  upon 
request,  on  a  test  basis,  to  all  Rescue  Coordination  Centers 
through  the  facilities  of  the  U.S.  Navy  Fleet  Numerical 
Weather  Facility,  Monterey,  California. 

2.  Input  Data 

a.   Leeway  Code 

For  this  program,  leeway  Is  expressed  as  a  two-digit 
code  as  shown  in  Table  1. 
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Leeway  Code 

00  No  Leeway 

10  Liferaft  with  Drogue 

11  Liferaft  without  Drogue 

12  Liferaft  (Minlmax) 

13  -  19   Downwind  drift  (minimax) 

Minimum  =  \%    of  wind  velocity 

Maximum  s  Second  digit  as  %   of 

> 
Wind  Velocity. 

20/30/etc.   Downwind  drift  at  2%,l%,  etc. 

of  the  wind  velocity 

21  -  99   Minimax  drift  at  some  fixed  percent 

of  wind  velocity  equal  to  first  digit; 

Deflection  to  the  right  and  left 

of  the  wind  line  equal  to 

10  degrees  times  the  second  digit. 

Table  1.   Leeway  Code 


Example:    kl    ■  k%   of  wind  velocity 
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+  30  degrees  off  wind  line 

b.  Number  of  Starting  Points 

The  number  of  starting  (or  distress)  positions  must 
be  indicated.  The  maximum  number  that  can  be  supplied  is 
n  i  ne. 

c.  Number  of  Datum  Times 

The  number  of  datum  (or  search)  times  for  which 
plots  are  requested  must  be  supplied.  The  maximum  number  of 
these  is  also  ni  ne. 

d.  Starting  Point  Data 

Starting  positions  must  be  entered  in  chronological 
order  and  each  point  must  be  f ul 1 v  identified  by: 

(1)  Date-time  group  (time  zone  Zulu) 

(2)  Month 

(3)  Latitude  and  Longitude  in  degrees  and  minutes. 
Tenths  of  minutes  must  be  specified.  The  number  of  these 
starting  positions  must  agree  with  paragraph  2. a. 

e.  Datum  Times 

Complete  Date-time  group  (including  month)  for  each 
datum  time  must  be  supplied/  in  chronological  order.  The 
number  of  these  datum  times  must  agree  with  paragraph  2.c. 

f.  Last  Known  Position 

This  is  the  last  known  position  of  the  search  object 
(see  Para.   611/  National   Search  and   Rescue  Manual)   to  be 
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used  in  computing  DR  position   error.   This  position  must  be 
complete  with  degrees,   minutes,  and  tenths  of   minutes.   If 
DR   position   error  is   not   applicable,   enter  one   of   the 
initial  positions  from  2.d. 
g.   Error  Factors 

These  factors  must  all  be  specified  to  at  least  one 
decimal  place,  in  the  following  order: 

(1)  Probable  navigation  error  of  search  craft. 

(2)  Initial  probable  error  in  position  of  the 
search  object. 

(3)  Percent  of   DR  distance   to  be   applied  as   DR 
error   (See  Para.  611,   National  Search  and    Rescue  Manual). 
Enter  0.0  if  not  applicable. 

h.   Search  Designators 

These  should  be  one-digit  numbers,  from  0  to  5, 
which  describe  which  Safety  Factors  are  to  be  applied  to  the 
Total  Probable  Error  of  position  for  each  datum  time.  One 
digit  must  be  supplied  for  each  datum  time. 

0  No  search  radius  is  to  be  computed  for 
that  datum  time 

1,2,3,4,5  The  Safety  Factor  corresponding  to 
search  1,2,3,4,  or  5  is  to  be  applied. 

3.   Message  Request  Format 

A  message  request   to  Fleet   Numerical  Weather   Central 
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will   insure  that   the  data   is  entered   correctly  into   the 
computer : 

FM  CEA 

TO  FLENUMWEACEN  MONTEREY 

BT 

UNCLAS 

DISTRESS  HOMEY  65  OVERDUE 

1.  REQUEST  SAR  PLOT. 

2.  DATA: 

A.  12/2/3// 

B.  030900Z  AUG/3520. 0N/05216.5W/031200Z  AUG, 
3600. 0N, 0511*0. 0W// 

C.  031600Z  AUG/032200Z  AUG/0U0200Z  AUG// 

D.  3  500.0N/051U0.0W/5.0/10.0/0.10/102// 

This  sample  message  requests  a  minimax  plot  of  a 
liferaft  from  two  initial  positions/  with  data  requested  for 
three  different  times  (para.  2. A.).  The  Initial  positions 
and  date-time  groups  are  given  in  para.  2.B.  (note  that  2 
digits  are  used  to  express  degrees  latitude  and  J5  digits  are 
used  for  degrees  longitude). 

Paragraph  2.C.  gives  the  date-time  groups  of  the  datum 
times.  Paragraph  2.D.  gives  the  last  known  position/  error 
factors/  and  search  designators.  for  this  example/  the 
first  search   is  to   be  computed   for  031600Z   and  a   second 
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search  for  0U0200Z.  No  search  radius  is  to  be  computed  for 
032200Z;  positions  and  position  errors  will  still  be 
prov  i  ded . 

k.      Caution:  Offshore  Regions 

This  program  is  designed  primarily  for  open  ocean 
search;  when  data  is  supplied  for  water  areas  within  100 
miles  from  shore/  answers  should  be  checked  by  hand 
calculations.  Requests  for  SAR  Plots  within  10  miles  from 
shore  should  be  avoided. 

Weather  information  used  in  this  program  is  applicable 
to  the  Northern  Hemisphere  only. 
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